Magnetic dilution in the geometrically frustrated SrCr 9p Gai 2 9p Oi9 and the role of 

local dynamics: a /iSR study 
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We investigate the spin dynamics of SrCr9 P Gai2-9 P Oig for 
p below and above the percolation threshold p c using muon 
spin relaxation. Our major findings are: (I) At T — > the 
relaxation rate is T independent and oc p 3 , (II) the slowing 
down of spin fluctuation is activated with an energy U which 
is also a linear function of p 3 and lim p _,o U = 8 K; this energy 
scale could stem only from a single ion anisotropy, and (III) 
the p dependence of the dynamical properties is identical be- 
low and above p c , indicating that they are controlled by local 
excitation. 

The SrCr 9p Gai2-9pOi9 [SCGO(p)] compound has 
been intensely studied experimentally in recent years, 
as it is believed to be a model compound for a 
"super-degenerate" antiferromagnct [jl] |],[l0| . By "super- 
degenerate" we mean that the classical ground state en- 
ergy is invariant under a local rotation of a small num- 
ber of spins, a symmetry which leads to local zero energy 
excitations in addition to the more common collective 
ones. One of the open questions in this area of research 
is which type of excitations will dominate in low tem- 
peratures. To date, this question was addressed only 
theoretically in two of the most famous super-degenerate 
magnets, the kagome and the pyrochlore where 
it was found that the dynamical properties are mostly 
controlled by the local excitations. Here we examine this 
question experimentally. We do so by measuring muon 
spin relaxation (/iSR) rates as a function of temperature 
T, magnetic field H, and, most importantly, magnetic ion 
concentration p above and below the percolation thresh- 
old p c . Our major finding is that the dynamical prop- 
erties of the system are impartial to p c , suggesting that 
they could not emerge from a collective phenomenon. 

Unexpectedly, our measurements also lead us to an- 
other finding regarding the spin Hamiltonian in SCGO. 
For long it was suspected that this Hamiltonian must 
contain terms other than the Heisenberg one. This is be- 
cause SCGO shows spin glass like effect in susceptibility 
experiments, such as a cusp at a temperature Tf ~ Ap K 
0, and a hysteresis between the field cooled (FC) and 
zero field cooled (ZFC) measurements. Neither of these 
could be understood in terms of Heisenberg spins on 



kagome or pyrochlore lattices [|13j . Our data indicate the 
presence of a single ion anisotropy with an energy scale 
of ~ 8 K which might be the origin of the susceptibility 
effects. 

Our samples were prepared by solid state reaction at 
1350C from the stoichiometric mixtures of Cr203, Ga203 
and SrC03. X-ray examination revealed the absence of 
foreign phase in the prepared samples and a slight smooth 
lattice expansion as p is decreased. Our intention was to 
have p values both above and below the p c of the lattice. 
However, there is a controversy whether SCGO repre- 
sents a kagome lattice or a pyrochlore-slab |3j . The p c for 
the kagome lattice is 0.6527 and for a pyrochlore-slab is 
not known. Nevertheless, the Curie- Weiss temperature 
8 CW in SCGO(p) is a linear function of p with a slope 
which is different below and above p = 0.61(5) in 
agreement with p c for the kagome lattice. Therefore we 
assume that p c — 0.6527 and prepare samples with p in 
the range 0.39 to 0.89. The value of p's in our samples is 
determined from their Curie constant. This method was 
found to be in agreement with the stoichiometric ratio 
in the sample preparation, the Curie- Weiss temperature, 
and Tf (when measurable) Q . 

Our fj.SH experiments were done in both TRIUMF and 
ISIS. In these experiments one follows the time evolu- 
tion of the spin polarization P z (t) of a muon implanted 
in a sample, through the asymmetry A(t) oc P z {t) in 
the positron emission of the muon decay. In addition, 
an external field H is applied along the initial muon 
spin (longitudinal) direction which defines the z axis. 
A(t) for three different samples at base temperature and 
H = 100 G is shown by the symbols in Fig. |l| where 
time is presented on a log scale. The small field of 100 G 
could be considered as zero field; it is applied in order to 
decouple the nuclear spin contribution to the muon spin 
relaxation. As can be seen from the figure, there is a 
strong variation in the time scale of relaxation between 
the different samples. In addition, the asymmetry for 
all samples has a "flat" beginning similar to Gaussian. 
In the inset of Fig. [l| we depict the asymmetry for the 
p = 0.89 at high temperature (5.5 K) where A is pre- 
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sented on a log scale. Clearly at high temperatures the 
relaxation is closer to exponential. Therefore, the asym- 
metry waveform is temperature dependent, as was first 
found in the p — 0.89 sample by Uemura et al. M, but 
depends weakly on the most important variable in this 
paper, namely, p. 

There are two ways to obtain a relaxation rate in a sit- 
uation where the waveform is changing: (I) using the 1/e 
criteria where we define the time T\ by A(T\) = A(Q)/e, 
or (II) fitting all data sets to a stretched exponential 

A(t)=A e X p(-(Xtf). (1) 

Representative fits are depicted in Fig. [| by the solid 
lines. However, in our case, where the stretched expo- 
nential fits the data well for more than 1/e of the initial 
asymmetry, as demonstrated by the arrow in Fig. [j], one 
finds that A = 1/T\. Therefore, both methods lead to 
the same relaxation rate although in (II) both A and (3 
are functions of T, iJ, and p whereas in (I) these param- 
eters impact only T\. We therefore continue the discus- 
sion using the stretched exponential fit approach which 
is more informative and accurate. Finally, in the samples 
with large p, Eq. |l] accounts only for early times where 
the asymmetry drops by ^80%. At later times a second 
component with low relaxation rate dominates. However, 
this second component is not observable in the small p 
samples due to experimental limitations imposed by the 
muon life time (2.2 fisec). Therefore, we concentrate here 
on the early time behavior of A(t). 

In Fig. U we show A in H = 100 G as a function of 
temperature for various values of p. All samples show 
critical slowing down starting at T — 20 — 5 K, followed 
by a saturation of the muon relaxation rate, in agree- 
ment with earlier /iSR work and more recent NMR 
measurements on the p — 0.89 sample. In Fig. ||a we 
depict (3 versus T for the different p's. As demonstrated 
before, the waveform in all our samples is exponential 
((3 — » 1) at high T but tends to be Gaussian ((3 — > 2) at 
low T. 

There are two possible mechanisms which can be re- 
sponsible for the loss of the muon polarization: static 
field distribution, or dynamic field fluctuations. It is 
possible to distinguish between these two by measuring 
the field dependence of the muon spin relaxation rate. 
In Fig. [| we show X(H ) on a semi-log scale for all our 
samples, and in Fig. ||b the field dependence of f3. The 
field dependence of A allows us to rule out the possibility 
of relaxation due to static field distribution as the fol- 
lowing argument shows. In the p = 0.89 sample, when 
H= 100 G, andT->0, the value of A is 10 /xsec -1 . If this 
relaxation would have been due to static field distribu- 
tion it would have implied an internal magnetic field [B] 
at the muon site in the order of 100 G (using A ~ 7 M [-B] 
where 7 M = 85.16 MHz/kG). The vector sum of this inter- 
nal field and an external longitudinal field of, say, 2 kG, 
would have been nearly parallel to the initial muon spin 
direction. Therefore, if the internal fields were static, we 



would expect a complete quenching of the relaxation rate 
in 2 kG and above, in contrast to observation. This line 
of argument applies to all other samples as well. Thus, 
the decay of P z (t) is not due to static random fields, and 
seems to be due to dynamic field fluctuations. 

On the other hand, the waveform is a Gaussian with 
a A that has a very weak field dependence for H up to 
2 kG; a field that obeys j^H/X 3> 1. This stands in 
strong contrast to all theories known to us of relaxation 
from dynamical fluctuations. These theories yield expo- 
nential relaxation when there is weak field dependence 
for j^H/X 3> 1. A dynamical Gaussian waveform is also 
very unusual experimentally, and is one of the on going 
puzzles of /iSR in SCGO @. Nevertheless, we interpret 
our data using a dynamical model, since the argument 
regarding the vector sum of internal and external fields 
seems more fundamental than the exact waveform. 

In dynamical models X(H) is proportional to the 
Fourier transform of the local field dynamical auto cor- 
relation function (Bi(t)Bj(0)) (i = x, y, z) at the fre- 
quency oj — 7 M -ff. We find that for all samples A (if) oc 
exp(— H/Hq) where iio = 1-16 T, as demonstrated by the 
solid line in Fig. [|. This fact indicates that the spectral 
density is not modified by magnetic dilution, apart from 
an over all factor, and that it is impartial to percolation. 

When no external field is applied the internal field fluc- 
tuation rate v is related to the muon relaxation rate by 

v<x{B 2 )X- x (2) 

where (-B 2 ) is the rms of the instantaneous field distri- 
bution at the muon site. We find that A _1 (T) could be 
well fitted to 

X- 1 (p,T) = Q(P) + Ccxpl-U( P )/T} (3) 

where C = 150(10) fisec is a global parameter. The fit 
results are given in Fig. || by the solid lines. This fit sug- 
gests that the internal field fluctuation rate is controlled 
by two dynamical processes: a temperature independent 
quantum one, and an activated one. A similar behav- 
ior was observed in the high spin molecul system CrNig 
where high T dynamics is due to over-the-barrier motion, 
and low T dynamics is due to quantum fluctuations [fL4| . 
Surprisingly, we find that Q~ l {p) and U{p) are linear 
functions of (p/p c ) 3 both below and above p c as demon- 
strated in Fig. ||. This result, together with the field 
dependence of A, leads us to our first main conclusion, 
namely, that the fluctuations are impartial to whether 
the lattice percolates and supports collective excitations 
or not. Therefore, the excitations are of local nature. 

It is interesting to compare this finding with other rele- 
vant experiments. Heat capacity measurements were per- 
formed by Ramirez et al. §. They found that C(T) - T 2 
even when the percolation threshold for the kagome lat- 
tice was crossed. They pointed out that this temper- 
ature dependence could result from collective antiferro- 
magnetic excitations in two dimensions of the acoustic 
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type, namely, u> oc k. A similar indication was made by 
Broholm et al. [[lj based on density of states p(ui) oc w 
found in neutron scattering. However, no dispersion re- 
lation of the acoustic type was ever found. Our finding 
of local excitations indicates that at low T the important 
w's are k independent. 

An extrapolation of Q^ 1 {p) in Fig. ||to p = suggests 
that Q diverges upon dilution. In fact for p < 0.05 we ex- 
pect Q(p) 3> C and the muon relaxation rate A should be 
T independent. In other words, SCGO(p < 0.05) should 
behave as if its spins are isolated. A similar extrapolation 
of U(p) to low p gives 8 K. This leads to our second major 
finding, namely, that the activated part of the dynamics 
does not emerge only from coupling between neighboring 
spins, but also from on site (single ion) interactions. The 
energy scale of this interaction is 8 K. When comparing 
this result with other experiments we find that it agrees 
with the energy scale of the anisotropy found by Schiffer 
et al. in their susceptibility measurements on SCGO sin- 
gle crystal ||, but disagrees with Ohta et al. who found 
a single ion energy two orders of magnitude smaller using 
EPR@. 

The conclusions drawn up to now are based on gross 
features in the data. Now we would like to speculate 
on how the p dependence of the muon relaxation rate 
A(p) is shared between the instantaneous field distribu- 
tion (-B 2 ) (p) and the fluctuation rate u(p) which deter- 
mines it using Eq. |^. First we calibrate (-B 2 ) (p) from 
the high temperature data where A shows a weak p de- 
pendence (See Fig. We assume that A cx p e for large 
T, where e is a small number. In addition, it is natural 
to expect v oc Jp 1 / 2 (where J is a coupling constant) in 
the high temperature range. Therefore, our calibration 
leads to {B 2 ) (p) oc p 1 / 2+£ , an expression which is only 
slightly different from the dilute limit where (-B 2 ) (p) oc p 
(e = 1/2). In SCGO the relation (B 2 ) (p) oc p 1 / 2+£ 
should hold at all temperatures since no static moment 
develops. Therefore, at base temperature, where A oc p 3 
we expect v oc p €_2 ' 5 , and since we can overrule small 
values of e there is a reasonable chance that v oc p~ 2 . 

In summary the dynamical spin fluctuations in SCGO 
are controlled by both quantum and activated dynamical 
processes. The activation energy is a linear function of p 3 , 
and indicates the presence of single ion anisotropy with 
an energy scale of 8 K. This anisotropy might be respon- 
sible for the spin glass like behavior of SCGO. The quan- 
tum fluctuation time (1/^) is most likely proportional to 
p . Both of the dynamical properties are impartial to the 
percolation threshold and indicate the local nature of the 
fluctuations. 
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FIG. 1. The asymmetry in the emission of muon decay 
positrons oc P z (t) as a function of time, for three different 
samples at base temperature, is presented on a semi-log scale 
by the symbols. The solid lines are fits to Eq. [l]. The in- 
set shows the asymmetry at high temperature, again on a 
semi-log scale. The solid line represents exponential decay. 

FIG. 2. Muon spin relaxation rates A as obtained from fits 
of raw data to Eq. |l| for various values of p. The solid lines 
are fits to Eq. ^. 

FIG. 3. The exponent /3 as obtained from a fit of the raw 
data to Eq. |l| for various values of p, (a) as a function of T at 
H — 100 G, and (b) as a function of H at base temperature. 

FIG. 4. The field dependence of the muon relaxation rate, 
plotted on a semi-log scale, for various values of p. The solid 
lines are parallel and are described in the text. 

FIG. 5. The parameters Q _1 and U of Eq. ^ plotted vs. 

{V/Pcf. 
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